****************************************************************************************
*** Paper: Populist Attitudes and Foreign Policy Postures
*** Tables A18-A21 in Appendix (Section "Robustness check: populism*nativism interaction")
*** Figures A3-A6 in Appendix  (Section "Robustness check: populism*nativism interaction")
****************************************************************************************

*** PROGRAM SETUP COMMANDS
	set more off            // Disable partitioned output
	set dp period			// Use comma instead of dot for decimals
	clear all               // Start with a clean slate
	set linesize 80         // Line size limit to make output more readable
	macro drop _all         // clear all macros
	set scheme s1mono, permanent			// set graph scheme, e.g., cleanplots, plottig, plotplainblind, burd3, s1mono
	graph set window fontface "Times New Roman"				// set font in graphs | "Arial" | "Times New Roman"

************************************************************************************************************************************
************************************************************************************************************************************

*** Open data set
	capture use "Data_clean.dta", clear


************************************************************************************************************************************

*** 1. Main effects models (Nativism instead of left-right)

	capture use "Data_clean.dta", clear


*** DV: Posture 1: Militant internationalism
	regress mi_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 1
	estimates store m1
	regress mi_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 2
	estimates store m2
	regress mi_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 3
	estimates store m3
	regress mi_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 4
	estimates store m4

*** DV: Posture 1: Cooperative internationalism
	regress ci_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 1
	estimates store m5
	regress ci_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 2
	estimates store m6
	regress ci_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 3
	estimates store m7
	regress ci_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 4
	estimates store m8


*** DV: Posture 1: Isolationism
	regress iso_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 1
	estimates store m9
	regress iso_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 2
	estimates store m10
	regress iso_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 3
	estimates store m11
	regress iso_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 4
	estimates store m12


*** DV: Posture 4: Global justice
	regress gj_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 1
	estimates store m13
	regress gj_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 2
	estimates store m14
	regress gj_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 3
	estimates store m15
	regress gj_index c.GoertzPopulism_0_1  nati_index age i.female i.educat2 if COUNTRY == 4
	estimates store m16


********************************************
*** Main effects tables (not reported)
********************************************
/*
	esttab 	m1 m2 m3 m4 using "Tables\Appendix_Table_A18.rtf", replace ///
		b(2) se(2) ar2 nogaps compress nobaselevels varwidth(15) ///
		modelwidth(6) nonumbers ///
		title("Regressions of militant internationalism posture on populism and other predictors") ///
		mtitles("France" "Germany" "UK" "Italy") ///
		drop()  ///
					coeflabels(_cons "Constant" ///
								GoertzPopulism_0_1 "Populism (Goertz)" ///
								nati_index "Nativism" ///
								age "Age" ///
								1.female "Gender (1=female)" ///
								2.educat2 "Education (2)" ///
								3.educat2 "Education (3)") ///
		order(_cons GoertzPopulism_0_1 nati_index age 1.female 2.educat2 3.educat2) ///
		note("Note: Reported are unstandardized OLS coefficients. All variables scaled 0-1.")


	esttab 	m5 m6 m7 m8 using "Tables\Appendix_Table_A19.rtf", replace ///
		b(2) se(2) ar2 nogaps compress nobaselevels varwidth(15) ///
		modelwidth(6) nonumbers ///
		title("Regressions of cooperative internationalism posture on populism and other predictors") ///
		mtitles("France" "Germany" "UK" "Italy") ///
		drop()  ///
					coeflabels(_cons "Constant" ///
								GoertzPopulism_0_1 "Populism (Goertz)" ///
								nati_index "Nativism" ///
								age "Age" ///
								1.female "Gender (1=female)" ///
								2.educat2 "Education (2)" ///
								3.educat2 "Education (3)") ///
		order(_cons GoertzPopulism_0_1 nati_index age 1.female 2.educat2 3.educat2) ///
		note("Note: Reported are unstandardized OLS coefficients. All variables scaled 0-1.")


	esttab 	m9 m10 m11 m12 using "Tables\Appendix_Table_A20.rtf", replace ///
		b(2) se(2) ar2 nogaps compress nobaselevels varwidth(15) ///
		modelwidth(6) nonumbers ///
		title("Regressions of isolationism posture on populism and other predictors") ///
		mtitles("France" "Germany" "UK" "Italy") ///
		drop()  ///
					coeflabels(_cons "Constant" ///
								GoertzPopulism_0_1 "Populism (Goertz)" ///
								nati_index "Nativism" ///
								age "Age" ///
								1.female "Gender (1=female)" ///
								2.educat2 "Education (2)" ///
								3.educat2 "Education (3)") ///
		order(_cons GoertzPopulism_0_1 nati_index age 1.female 2.educat2 3.educat2) ///
		note("Note: Reported are unstandardized OLS coefficients. All variables scaled 0-1.")


	esttab 	m13 m14 m15 m16 using "Tables\Appendix_Table_A21.rtf", replace ///
		b(2) se(2) ar2 nogaps compress nobaselevels varwidth(15) ///
		modelwidth(6) nonumbers ///
		title("Regressions of global justice posture on populism and other predictors") ///
		mtitles("France" "Germany" "UK" "Italy") ///
		drop()  ///
					coeflabels(_cons "Constant" ///
								GoertzPopulism_0_1 "Populism (Goertz)" ///
								nati_index "Nativism" ///
								age "Age" ///
								1.female "Gender (1=female)" ///
								2.educat2 "Education (2)" ///
								3.educat2 "Education (3)") ///
		order(_cons GoertzPopulism_0_1 nati_index age 1.female 2.educat2 3.educat2) ///
		note("Note: Reported are unstandardized OLS coefficients. All variables scaled 0-1.")
*/


************************************************************************************************************************************

************************************************************************************************************************************

*** 2. Interaction models

	capture use "Data_clean.dta", clear


*** DV: Posture 1: Militant internationalism
	regress mi_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 1
	estimates store m1
	regress mi_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 2
	estimates store m2
	regress mi_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 3
	estimates store m3
	regress mi_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 4
	estimates store m4

*** DV: Posture 1: Cooperative internationalism
	regress ci_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 1
	estimates store m5
	regress ci_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 2
	estimates store m6
	regress ci_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 3
	estimates store m7
	regress ci_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 4
	estimates store m8


*** DV: Posture 1: Isolationism
	regress iso_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 1
	estimates store m9
	regress iso_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 2
	estimates store m10
	regress iso_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 3
	estimates store m11
	regress iso_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 4
	estimates store m12


*** DV: Posture 4: Global justice
	regress gj_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 1
	estimates store m13
	regress gj_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 2
	estimates store m14
	regress gj_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 3
	estimates store m15
	regress gj_index c.GoertzPopulism_0_1##c.nati_index  age i.female i.educat1 if COUNTRY == 4
	estimates store m16


**********************************
*** Tables A18-A21 in Appendix
**********************************

	esttab 	m1 m2 m3 m4 using "Tables\Appendix_Table_A18.rtf", replace ///
		b(2) se(2) ar2 nogaps compress nobaselevels varwidth(15) ///
		modelwidth(6) nonumbers ///
		title("Regressions of militant internationalism posture on populism*nativism and other predictors") ///
		mtitles("France" "Germany" "UK" "Italy") ///
		drop()  ///
					coeflabels(_cons "Constant" ///
								GoertzPopulism_0_1 "Populism (Goertz)" ///
								nati_index "Nativism" ///
								c.GoertzPopulism_0_1#c.nati_index "Populism * Nativism" ///
								age "Age" ///
								1.female "Gender (1=female)" ///
								2.educat1 "Education (2)" ///
								3.educat1 "Education (3)") ///
		order(_cons GoertzPopulism_0_1 nati_index c.GoertzPopulism_0_1#c.nati_index age 1.female 2.educat1 3.educat1) ///
		note("Note: Reported are unstandardized OLS coefficients. All variables scaled 0-1.")

	esttab 	m5 m6 m7 m8 using "Tables\Appendix_Table_A18.rtf", replace ///
		b(2) se(2) ar2 nogaps compress nobaselevels varwidth(15) ///
		modelwidth(6) nonumbers ///
		title("Regressions of cooperative internationalism on populism*nativism and other predictors") ///
		mtitles("France" "Germany" "UK" "Italy") ///
		drop()  ///
					coeflabels(_cons "Constant" ///
								GoertzPopulism_0_1 "Populism (Goertz)" ///
								nati_index "Nativism" ///
								c.GoertzPopulism_0_1#c.nati_index "Populism * Nativism" ///
								age "Age" ///
								1.female "Gender (1=female)" ///
								2.educat1 "Education (2)" ///
								3.educat1 "Education (3)") ///
		order(_cons GoertzPopulism_0_1 nati_index c.GoertzPopulism_0_1#c.nati_index age 1.female 2.educat1 3.educat1) ///
		note("Note: Reported are unstandardized OLS coefficients. All variables scaled 0-1.")

	esttab 	m9 m10 m11 m12 using "Tables\Appendix_Table_A18.rtf", replace ///
		b(2) se(2) ar2 nogaps compress nobaselevels varwidth(15) ///
		modelwidth(6) nonumbers ///
		title("Regressions of isolationism posture on populism*nativism and other predictors") ///
		mtitles("France" "Germany" "UK" "Italy") ///
		drop()  ///
					coeflabels(_cons "Constant" ///
								GoertzPopulism_0_1 "Populism (Goertz)" ///
								nati_index "Nativism" ///
								c.GoertzPopulism_0_1#c.nati_index "Populism * Nativism" ///
								age "Age" ///
								1.female "Gender (1=female)" ///
								2.educat1 "Education (2)" ///
								3.educat1 "Education (3)") ///
		order(_cons GoertzPopulism_0_1 nati_index c.GoertzPopulism_0_1#c.nati_index age 1.female 2.educat1 3.educat1) ///
		note("Note: Reported are unstandardized OLS coefficients. All variables scaled 0-1.")

	esttab 	m13 m14 m15 m16 using "Tables\Appendix_Table_A18.rtf", replace ///
		b(2) se(2) ar2 nogaps compress nobaselevels varwidth(15) ///
		modelwidth(6) nonumbers ///
		title("Regressions of global justice posture on populism*nativism and other predictors") ///
		mtitles("France" "Germany" "UK" "Italy") ///
		drop()  ///
					coeflabels(_cons "Constant" ///
								GoertzPopulism_0_1 "Populism (Goertz)" ///
								nati_index "Nativism" ///
								c.GoertzPopulism_0_1#c.nati_index "Populism * Nativism" ///
								age "Age" ///
								1.female "Gender (1=female)" ///
								2.educat1 "Education (2)" ///
								3.educat1 "Education (3)") ///
		order(_cons GoertzPopulism_0_1 nati_index c.GoertzPopulism_0_1#c.nati_index age 1.female 2.educat1 3.educat1) ///
		note("Note: Reported are unstandardized OLS coefficients. All variables scaled 0-1.")


*** Figures militant internationalism
	estimates restore m1
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g1a, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Militant internationalism", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("France", size(large)) ///
								legend(off) ///
								text(.85 -0.01 "Nativism max", color(gs9) place(right) size(large)) ///
								text(.40 -0.01 "Nativism min", color(gs0) place(right) size(large))

	estimates restore m2
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g1b, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Militant internationalism", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("Germany", size(large)) ///
								legend(off) 

	estimates restore m3
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g1c, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Militant internationalism", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("UK", size(large)) ///
								legend(off)

	estimates restore m4
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g1d, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Militant internationalism", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("Italy", size(large)) ///
								legend(off) 


*** Figures Cooperative internationalism
	estimates restore m5
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g2a, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Cooperative internationalism", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("France", size(large)) ///
								legend(off) ///
								text(.30 -0.01 "Nativism max", color(gs9) place(right) size(large)) ///
								text(.85 -0.01 "Nativism min", color(gs0) place(right) size(large))

	estimates restore m6
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g2b, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Cooperative internationalism", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("Germany", size(large)) ///
								legend(off) 

	estimates restore m7
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g2c, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Cooperative internationalism", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("UK", size(large)) ///
								legend(off)

	estimates restore m8
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g2d, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Cooperative internationalism", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("Italy", size(large)) ///
								legend(off) 


*** Figures Isolationism
	estimates restore m9
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g3a, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Isolationism", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("France", size(large)) ///
								legend(off) ///
								text(.80 -0.01 "Nativism max", color(gs9) place(right) size(large)) ///
								text(.10 -0.01 "Nativism min", color(gs0) place(right) size(large))

	estimates restore m10
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g3b, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Isolationism", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("Germany", size(large)) ///
								legend(off) 

	estimates restore m11
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g3c, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Isolationism", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("UK", size(large)) ///
								legend(off)

	estimates restore m12
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g3d, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Isolationism", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("Italy", size(large)) ///
								legend(off) 



*** Figures Global justice
	estimates restore m13
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g4a, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Global justice", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("France", size(large)) ///
								legend(off) ///
								text(.10 -0.01 "Nativism max", color(gs9) place(right) size(large)) ///
								text(.80 -0.01 "Nativism min", color(gs0) place(right) size(large))

	estimates restore m14
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g4b, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Global justice", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("Germany", size(large)) ///
								legend(off) 

	estimates restore m15
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g4c, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Global justice", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("UK", size(large)) ///
								legend(off)

	estimates restore m16
						margins, at(GoertzPopulism_0_1=(0 (.025) 1) nati_index=(0 1))
						marginsplot,  ///
								name(g4d, replace) ///
								recast(scatter) ///
								plotopts(mcolor(none))  ///
								ciopts(recast(rbar) barwidth(0.01)) ci1opts(color(gs0) alcolor(none)) ci2opts(color(gs9%80) alcolor(none)) ///
								xtitle("Populism (Goertz)", size(large)) xscale(range(0 1)) xlabel(0 (0.2) 1, labsize(large)) ///	
								ytitle("Global justice", size(large)) yscale(range(0 1)) ylabel(0 (.2) 1, labsize(large)) ///	
								title("Italy", size(large)) ///
								legend(off) 
